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Abstract 

Semi-empirical method of calculation of quenching factors for scintillators is de- 
scribed. It is based on classical Birks formula with the total stopping powers for 
electrons and ions which are calculated with the ESTAR and SRIM codes, respec- 
tively. Method has only one fitting parameter (the Birks factor kB) which can have 
different values for the same material in different conditions of measurements and 
data treatment. A hypothesis is used that, once the kB value is obtained by fitting 
data for particles of one kind and in some energy region (e.g. for a few MeV a 
particles from internal contamination of a detector), it can be applied to calculate 
quenching factors for particles of another kind and for another energies (e.g. for 
low energy nuclear recoils) if all data are measured in the same experimental condi- 
tions and are treated in the same way. Applicability of the method is demonstrated 
on many examples including materials with different mechanisms of scintillation: 
organic scintillators (solid CsHg, and liquid CigHis, C9H12); crystal scintillators 
(pure CdW0 4 , PbW0 4 , ZnW0 4 , CaW0 4 , CeF 3 , and doped CaF 2 (Eu), CsI(Tl), 
CsI(Na), Nal(Tl)); liquid noble gases (LXe). Estimations of quenching factors for 
nuclear recoils are also given for some scintillators where experimental data are 
absent (CdW0 4 , PbW0 4 , CeF 3 , Bi 4 Ge 3 0i 2 , LiF, ZnSe). 
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1 Introduction 



In accordance with our current understanding of astronomical observations, 
usual matter constitutes only ~ 4% of the Universe; the main components 
are dark matter (~ 23%) and dark energy (~ 73%) [lj. Various extensions 
of the Standard Model propose many candidates on the role of dark matter 
(DM) particles [2] which are neutral and only weakly interact with matter 
(Weakly Interacting Massive Particles, WIMPs). One of the approaches to 
discover these particles is to detect scattering of WIMPs on atomic nuclei 
in sensitive detectors placed deep underground and measured in extra low 
background conditions [3]. Taking into account likely mass range and velocities 
of WIMPs, energies of nuclear recoils are expected below ~ 100 keV with 
character interaction rates of 1 — 10~ 6 events kg -1 d _1 . Many searches of 
WIMPs with semiconductor, scintillator and bolometer detectors to-date gave 
only negative results (see [3] and references therein); instead positive evidence 
for DM particles (WIMPs are a subclass; other candidates and other kinds of 
interactions are also available) in the galactic halo has been pointed out by 
DAMA experiments by exploiting the DM annual modulation signature with 
Nal(Tl) scintillators during more than 10 years long measurements [I]. 

For a long time it is known that amount of light produced in scintillating 
material by highly ionizing particles is lower than that produced by electrons 
of the same energy [5]. Thus, in a scintillator calibrated with electron and/or 
7 sources (which is an usual practice), signals from ions will be seen at lower 
energies (sometimes up to ~ 40 times) than their real values. Evidently knowl- 
edge of these transformation coefficients - quenching factors - is extremely 
important in searches for WIMPs and in predictions where the WIMPs signal 
should be expected. Many experimental efforts were devoted to measurements, 
sometimes very sophisticated, of quenching factors at low energies in different 
detectors (see e.g. recent works pi7TI8f5][TUjlllll2] and further references). 

Quenching factors are also needed in measurements and interpretation of sig- 
nals from a particles in scintillators. As examples, we can mention here recent 
experiments on searches (and first observations) of extremely rare a decays 
(Ti /a = 10 18 - 10 19 yr): 180 W in CdW0 4 [13] and CaW0 4 [H crystal scintil- 
lators, and 151 Eu in CaF 2 (Eu) [15]. 

While few approaches in calculation of quenching factors are known [i~6fl7|ll8|[T9l 
satisfactory theory able to exactly predict (and very often even to describe 
already measured) quenching factors for all detectors and particles still is ab- 
sent. For example, in the Lindhard's approach [18p1 it is possible to calculate 
quenching of ions with atomic number Z in scintillator only with the same Z 



2 Good description is given in more accessible source [20 
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number; in addition, this theory predicts decrease of quenching factors at low- 
energies, very often in contradiction with experimental data. Hitachi's model 
[T9] gives better description and for wider data range, however it is not easy 
to reproduce these calculations independently. 

Below we describe rather simple method of calculation of quenching factors 
for different ions (from protons to heavy recoils), based on semi-empirical 
approach of Birks [16] and using available in Internet software for calcula- 
tion of stopping powers for electrons and ions (ESTAR [21J and SRIM [22] 
codes, respectively). It employs only one parameter (kB Birks factor) which 
could be found by fitting experimental data measured for particles of one kind 
in some energy region (e.g. for a particles from external sources or internal 
contamination of a detector by U/Th chains, 147 Sm, 190 Pt, etc.) but after- 
wards can be used to calculate quenching factors for other particles and in 
other energy regions (e.g. for nuclear recoils at low energies). Summary of the 
method is given in section 2. Calculations with this method are demonstrated 
in section 3 for number of scintillators: organic scintillators (solid CsHg, and 
liquid Ci 6 Hi 8 , C 9 Hi 2 ); crystal scintillators (pure CdW0 4 , PbW0 4 , ZnW0 4 , 
CaWC-4, CeF 3 , and doped CaF 2 (Eu), CsI(Tl), CsI(Na), Nal(Tl)); liquid noble 
gases (LXe). Estimations of quenching factors for nuclear recoils are also given 
for some scintillators where experimental data are absent (CdW04, PbW04, 
CeF 3 , Bi 4 Ge30i2, LiF, ZnSe). Section 4 gives conclusions. 



2 Outlines of the method 

In calculation of quenching factors, we follow Birks approach in description 
of quenching of the light yield for highly ionizing particles PfT6] . Light yield 
of scintillating material depends not only on energy of particle E but also on 
how big is its stopping power dE/dr in the material. Fig. 1 gives example of 
stopping powers for electrons (calculated with the ESTAR software [2T]). and 
for protons, alpha particles, O, Ca and W ions (calculated with the SRIM 
code [22]) in the CaW04 material. 

In case when created in a scintillator excitation centers are spaced at large 
distances and interactions between them can be neglected, what is realized 
for particles with low stopping power (fast electrons, energies above E ~ 100 
keV), scintillation yield dL is proportional to released energy dE: dL = SdE 
(where S is the absolute scintillation factor), or in differential form 

dr dr 

To account for suppression of the light yield for highly ionizing particles (pro- 
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Fig. 1. (Color online) Stopping powers for electrons, protons, a particles, O, Ca and 
W ions in CaW04. In (a) only total dE/dr are given; in (b) also nuclear (dotted 
line) and electronic (dashed line) parts of dE/dr are drawn for a particles and W 
ions; for electrons, collision and radiation parts are shown. 

tons, a particles and nuclear recoils; hereafter all of them will be named 
"ions"), Birks proposed semi-empirical formula [5fl~6] : 



where BdE/dr is density of excitation centers along the track, and k is a 
quenching factor; kB is usually treated as a single parameter (Birks factor). 

Equation (2) gives the following approximations for light yields for particles 
with low (fast electrons) and high (ions) stopping power: 



dL 

dr 
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but in general light yield is: 




(4) 



Quenching factor for ion is a ratio of light yield of ions to that of electrons 



3 Quenching factor for a particles is often named "a//5" ratio. 



of the same energy: 



ft(^) = ^ = 4™. (5) 



Le(E) 



The 5 factor disappeared in the ratic0, and Qi(E) depends only on single 
parameter kB. 

Instead of quenching factor, sometimes a relative light yield, ratio of Li to 
energy E, normalized to that for electron L e at some energy E , is used: 



Relation between Qi(E) and Ri(E) is evident: 

R(F] _ HE) L e (E)/E _ L e (E)/E 



Thus, Ri(E) is close to Qi(E) if E and Eq are in energy region where electron 
light yield L e is proportional to energy. 

Taking into account that approximately (see Eq. (3)) dL e /dE = S and dLi/dE - 
W (dE^dr) ' we can °btain the following approximation for quenching factor: 

. fm U{E) L t (E)/E dU/dE 1 
Qi\E) 



L e {E) L e {E)/E dL e /dE kB(dE/dr)i' v ' 

While this expression^ is approximate, it gives the following important fea- 
tures of quenching factor: 

(1) Quenching factor depends on energy. This is not so trivial feature because 
in many papers on Qi measurements it was supposed that Qi is constant; 

(2) Qi is minimal when (dE/dr)i is maximal; 

(3) Qi increases at low energies; this is a consequence of decrease of (dE/dr)i, 
see Fig. 1. 



Thus, the S factor is supposed independent on energy and equal for electrons 
and ions; in the following we will suppose that the kB factor is also independent on 
energy. 

5 In [5j it was erroneously written as Qi = {dE / dr)i/kB, see Eq. (6.5). 
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In the following, we will use Eq. (5) to calculate quenching factors for different 
particles and different scintillators. Stopping powers for ions will be calculated 
with the SRIM code [22] • It should be noted that it is possible to calculate 
stopping powers for a particles and protons also with the ASTAR and PSTAR 
codes of the STAR package |21] , respectively, but list of materials available 
is restricted. There is also difference in dE/dr calculated with the SRIM and 
ASTAR & PSTAR for a particles and protons which results also in difference 
in calculated quenching factors; some examples are given later. The SRIM 
code does not allow to calculate stopping powers for electrons, and for this 
the ESTAR code [21] will be used. Contrary to other approaches, we will use 
total stopping powers instead of using only electronic part of dE/dr. Currently 
we accept this as a hypothesis and will show in future that it works well (in 
particular, see Fig. 9 later). 

Before to calculate quenching factors and compare them with experimental 
values measured in different works, the following general note should be made. 
Quenching factors could depend on many conditions of experimental measure- 
ments: 

(1) If scintillator is not pure but doped with some material which enhances 
its scintillating characteristics (e.g. Tl in Nal(Tl) or CsI(Tl)), Qi depends on 
kind and amount of dopant. For example, for PbW04 detectors and external 
a particles of 5.25 MeV Qi were measured as Qi = 0.19 — 0.32 with different 
dopants and their different amounts [23]- And even if some material is con- 
sidered as a "pure" scintillator, usually it also contains impurities and defects 
which could affect Qi. 

(2) Scintillation dynamics and light output depend on temperature. This is 
also well known experimental fact; f.e. we can quote Ref. [21] where average 
decay time of CaMo04 scintillator was measured as ~ 17 /is at +20° C but as ~ 
350 /is at —140° C. Usually scintillation has few components with amplitudes 
different for different particles. Temperature dependence of these amplitudes 
could lead to change in quenching factors; for example, for a particle with 
E a = 2.14 MeV emitted in a decay of 152 Gd inside Gd 2 Si0 5 (Ce) detector, 
the a/P ratio changed from 0.168 at -20° C to 0.178 at +20° C [25]. In Ref. 
[2S] , change in temperature from +20° C to —20° C resulted in increase of the 
a/p ratio on 7% in Nal(Tl), 35% in CsI(Tl) and 25% in CsI(Na); change from 
+20° C to +80° C decreased a/P ratio on 3%, 15% and 30% in these crystals, 
respectively. 

(3) Such a technical parameter as time At during which scintillation signal 
is collected by a data acquisition system, is in fact very important and could 
drastically change Qi values. This is because the amplitudes of different com- 
ponents of scintillation signal are different for different particles; thus different 
parts of a signal will be collected during At. This is illustrated by Fig. 2 where 
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relative light output for CsI(Tl) and external a particles with energies between 
~ and 10 MeV were measured with At = 1 /is and 7 /is |27J. Attempts to fit 
these data also are shown. Descriptions of the data with Eq. (6) (which is bet- 
ter for At = 7 /is) lead to different values of the kB parameter: kB = 1.1 x 10~ 3 
g MeV -1 cm" 2 for At = 1 /js, and kB = 2.3 x 1CT 3 g MeV" 1 cm" 2 for At = 7 
/j,s. This gives very important conclusion: the Birks factor kB, which very often 
is named "Birks constant" , in fact is not a fundamental constant for a given 
material but could have different values at different experimental conditions 
(including time of a signal collection At). 




2000 4000 6000 8000 10000 

Energy (keV) 

Fig. 2. (Color online) Relative light output (normalized to that for electrons of 662 
keV) for CsI(Tl) and a particles [27J. Fits of the data with Eq. (6) are shown as 
continuous curves calculated with the ESTAR (for electrons), and ASTAR or SRIM 
(for a's) codes. 

If signals are not collected during a proper time, it is possible to obtain wrong 
conclusions on Qi values. For example, the light output for protons of 662 keV 
in CsI(Tl) measured during At = 1 /is in [27] is higher than that for 7 quanta 
of the same energy because scintillation signal for p is faster. Thus, instead of 
expected quenching (Qi < 1) we, on contrary, obtain enhancement (Qi > 1). 
However, with 7 and p signals collected during longer time of 7 /xs, we have 
the usual situation when Qi < 1 [27] . 

(4) Sometimes quenching factors are derived from measurements with non- 
monoenergetic neutron sources like Am-Be[£]- After collision with neutron, 
nuclear recoil with mass M r has energy determined by initial energy of neutron 
E n and neutron scattering angle 9: 

2E 1 2E 

E r = ^ /i + s in 2 #-cosfl v V -sin 2 9) ~ — (1 -cos 6), (9) 

1 + itr a 



6 Spectrum of neutrons from Am-Be (Am-B) source has complex structure with 
energies up to ~ 11 MeV (~ 5.5 MeV) [28]. 
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where /i = M r /m n , m n is the neutron mass, and the last approximation is 
valid for heavy nuclei (// ^> 1). Energy dependence of quenching factor can 
be found only if the initial neutron energy and scattering angle are known; 
measurements with non-monoenergetic neutron sources give only some average 
Qi value which, nevertheless, could be useful estimate of Qi for energy range 
~ — 4E n /A, where A is mass number of nuclear recoil, if more detailed data 
are absent (this effective energy range is also related with dependence of the 
cross section on angle of scattering). 

(5) Presence of electric field in case of liquid noble gases could distort initial 
Qi values (obtained without electric field). Other phenomena like channeling 
effect in crystals (see f.e. [29]), dependence of Qi on direction of particle's 
movement relatively to crystal axes (as f.e. in CdW04 [13]) or diffusion and 
movement of molecules in liquids [SU] also effect quenching factors. 

In accordance with the above mentioned, we will not expect that the Birks 
factor kB for a given material will have the same value in different measure- 
ments. However, we will expect that if conditions of measurements and data 
treatment are fixed, kB will be the same for all particles. Such a hypothesis 
was discussed already in [5] , and it was supported by some experimental data. 
Below we show that it gives reliable results for a range of energies of interest 
here (low energy ions and a particles with energies up to ~ 10 MeV). 



3 Calculation of quenching factors 

Results of calculation of quenching factors with Eq. (5), or relative light out- 
puts with Eq. (6) are presented below for a number of organic, crystal and 
liquid noble gases scintillators. To compare calculations with experimental 
results, among big number of experimental papers we mainly chose more re- 
cent articles where data were obtained with better techniques (monoenergetic 
neutron beam instead of Am-Be source, for example) or/and in wider energy 
range. 

3.1 Organic scintillators 
3.1.1 Polystyrene (C S H 8 ) 

Relative light output L a (E)/E for a particles with energies of 2 — 9 MeV 
(normalized to that for electrons of 976 keV0) for polystyrene scintillator 

7 Difference between the relative light output and quenching factor for presented 
energies is small. 
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(chemical formula CsHg, density p = 1.06 g cm~ 3 ) was measured recently in 
Ref. [31] using external a particles from 241 Am source with a set of thin mylar 
absorbers, and internal a particles from contamination of the scintillator by 
U/Th chains. Fit of these experimental data by calculations with Eq. (6) and 
with the Birks factor kB = 9.0 x 10~ 3 g MeV -1 cm 2 is presented in Fig. 3a. 
Fitting curves were obtained with stopping powers calculated with the SRIM 
or ASTAR codes for a particles, and ESTAR code for electrons. Difference 
between the ASTAR and SRIM calculations is not very big, however x 2 /n.d.f. 
vah43 for the ASTAR (0.50) is better than that for the SRIM (0.58). 

3.1.2 PXE (C 16 H W ) 

In Ref. [32] , quenching factors for a particles were investigated for liquid scin- 
tillator: phenyl-o-xylylethane (l,2-dimethyl-4-(l-phenylethyl)-benzene, PXE, 
chemical formula CigHig, density p = 0.988 g cm -3 ) doped with para-Terphenyl 
(1,4-diphenylbenzene, p-Tp) at 2.0 g/1 and bis-MSB (l,4-bis(2-methylstyryl)- 
benzene) at 20 mg/1. Alpha particles belong to internal contamination of 
scintillator by daughters from 238 U chain ( 222 Rn, 218 Po, 214 Po, 210 Po); E a = 

5.3 - 7.7 MeV. 

In fact, in case of a decay of nucleus inside a scintillator, released light output 
has two components: from a particle and from nuclear recoil. However, cor- 
rections for the light output from heavy nuclear recoils are small: for example, 
for a decay of 210 Po (Q a = 5407 keV) corresponding energies are E a = 5304 
keV, E r = 103 keV and with the value of kB = 6.8 x 10~ 3 g MeV" 1 cm" 2 (see 
Fig. 3b) the light output from 206 Pb recoil is only 2.1% of that from the a 
particle. We will neglect these corrections in the following. Quenching factors 
calculated for Ci 6 Hi 8 (PXE) with kB = 6.8 x 10 -3 g MeV -1 cm -2 are shown 
in Fig. 3b in comparison with the experimental data |32j (% 2 /n.d.f. = 0.08). 

3.1.3 Pseudocumene (C 9 H 12 ) 

(1) In Ref. [32], quenching factors for a particles were studied as above also for 
another liquid scintillator: pseudocumene (1,2,4-trimethylbenzene, PC, C9H12, 
p = 0.876 g cm -3 ) doped with 2,5-diphenyloxazole (PPO) at a concentration 
of 1.5 g/1. Fit of these data (x 2 /n.d.f. = 0.12) is shown in Fig. 3c with kB = 

9.4 x 10~ 3 g MeV- 1 cm" 2 . 

(2) Relative light outputs Li(E)/E for protons and C ions in pseudocumene 
(BC505 liquid scintillator) were measured in [33] in range of energies of 29—943 
keV (protons) and 46 — 500 keV (C ions). They were normalized to relative 

8 Values of x 2 / n -d.f. are calculated everywhere without taking into account uncer- 
tainties in energy. 
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Fig. 3. (Color online) Quenching factors or relative light outputs for a particles and 
ions in organic scintillators: experimental data and their fit with Eq. (5) or (6). (a) 
C 8 H 8 scintillator (polystyrene) [HI], kB = 9.0 x 1(T 3 g MeV" 1 cm" 2 ; (b) Ci 6 Hi 8 
liquid scintillator (PXE) [32], kB = 6.8 x 1(T 3 g MeV" 1 cm' 2 ; (c) C 9 Hi 2 liquid 
scintillator (pseudocumene) [32], kB = 9.4 x 1CT 3 g MeV" 1 cm" 2 ; (d) C 9 H 12 [33], 
kB = 35 x 1CT 3 g MeV" 1 cm" 2 ; E = 22 keV for protons and 32 keV for C ions. In 
part (a), results with stopping powers calculated with the ASTAR and SRIM codes 
for a particles are shown. 

light output for electrons at 22 keV (for p) and 32 keV (for C) [33]. The Birks 
factor kB = 42 x 10 -3 g MeV -1 cm -2 was found by fitting the experimental 
data for protons by Eq. (6) (see Fig. 3d); xV n -d.f. value is 2.8, with main 
contribution due to the last experimental point. After this, the curve for C 
ions was calculated with this kB value. Comparison of the curve with the C 
data gives x 2 / n -d.f. = 5.3. While this value is high, nevertheless calculations 
are in a proper agreement with the measured data, as one can see in Fig. 3d. 



The kB values for pseudocumene C 9 H 12 obtained by fitting data of |32j and 
[33] . respectively (Fig. 3c and 3d), are quite different. However, this is not a 
surprise taking into account different conditions of measurements and probably 
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dopants used. Nevertheless, the same kB value fixed in one experiment [33] 
allowed to describe both data sets: for protons and C ions. 

3.2 Crystal scintillators 

3.2.1 CdWO A 

(1) Quenching factors for a particles in CdWC>4 crystal scintillator (density 
p = 7.9 g cm -3 ) were measured in experimental searches for rare a decay of 
180 W (Q a = 2516 keV) with CdW0 4 detector [T3] where it was observed at the 
first time (Ty 2 — 1.1 x 10 18 yr). Alpha particles from external 241 Am source 
(with a set of thin absorbers 9 1) and from internal contamination of CdWC>4 
by U/Th chains were used; E a = 0.47 — 8.79 MeV. The experimental points 
and fitting curve calculated with Eq. (5) and kB = 10.1 x 10~ 3 g MeV -1 cm -2 
are shown in Fig. 4a (x 2 /n.d.f. = 13; without last point xV n -d.f. = 2.2). 




a particles 



2000 4000 6000 8000 10000 1000 2000 3000 4000 5000 6000 

Energy (keV) Energy (keV) 

Fig. 4. (Color online) Quenching factors for a particles and protons in CdW04. (a) 
Experimental data from Ref. [13] and their fit by Eq. (5) with kB = 10.1 x 10~ 3 
g MeV" 1 cm -2 ; (b) data for protons and a's from [34] and their fit with 
kB = 21.5 x 10~ 3 g MeV" 1 cm" 2 . 

(2) In Ref. [31], quenching factors for protons provided by accelerator with 
E p = 2.8 MeV and 3.4 MeV were measured as 0.34 ± 0.02 and 0.39 ± 0.02, 
respectively. Fit of these points by Eq. (5) (x 2 /n.d.f. = 0.90) was possible with 
the value of the Birks factor: kB = 21.5 x 10~ 3 g MeV -1 cm -2 , very different 
from that obtained for CdW04 by fitting data from Ref. [13] . Quenching fac- 
tors for a particles were also measured in [34J. A mixed source with a's from 
239 Pu (E a ~ 5.1 MeV [35]), 241 Am (E a ~ 5.5 MeV) and 244 Cm (E a ~ 5.8 

9 It should be noted that quenching factors in CdW04 depend on direction of 
movement of a particle; we use here data for direction perpendicular to (010) crystal 
plane measured in wider energy range [T3] . 
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MeV) was used. Quenching factor of Q a ~ 0.1 was measured (different from 
those in [13J). The value of kB = 21.5 x 10~ 3 g MeV -1 cm -2 obtained for pro- 
tons gives theoretical curve for a's in good agreement with this experimental 
result (x 2 /n.d.f. = 0.12), see Fig. 4b. 

3.2.2 CaF 2 (Eu) 

(1) Quenching factors for a particles in CaF2 crystal scintillator (p = 3.18 g 
cm -3 ) doped by Eu at 0.4% were measured in [15] with external ( 241 Am) and 
internal (U/Th chains, 147 Sm) a sources (E a = 1 — 9 MeV). The measured 
experimental points and their fit by Eq. (5) with the Birks factor kB = 5.3 x 
10~ 3 g MeV" 1 cm -2 are presented in Fig. 5a. Stopping powers of a particles 
were calculated with the SRIM and ASTAR codes; fitting with the ASTAR 
curve is better: % 2 /n.d.f. is 3.5 (with near 50% contribution from the last point) 
while for the SRIM % 2 /n.d.f. = 7.6. It is a pity that the STAR package allows 
calculations of dE/dr only for restricted list of materials in case of a particles 
and protons (and other ions are absent). From the other side, SRIM allows to 
calculate dE/dr for any materials and for any ions, but not for electrons. Fig. 
5a gives an idea that probably calculations of dE/dr for all particles inside 
the same package would give better description of quenching factors. 




a particles 
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Fig. 5. (Color online) Quenching factors for a particles, and F and Ca ions in 
CaF2(Eu). (a) Experimental data from Ref. [15] for a particles and their fit by Eq. 
(5) with kB = 5.3 x 10~ 3 g MeV -1 cm -2 with stopping powers calculated with the 
SRIM and ASTAR. (b) Data for F (squares) and Ca (triangles) ions from |36| and 
their fit with kB = 5.1 x 10" 3 g MeV" 1 cm . (c) Data from [37J for a's with the 
same kB value as in (b). 

(2) Quenching factors for F and Ca ions in CaF 2 (Eu) measured in [3B] are 
shown in Fig. 5b. The kB value was obtained by fitting the F data as: kB = 
5.1 x 10~ 3 g MeV -1 cm" 2 (% 2 /n.d.f. = 1.3). Curve for Ca ions was calculated 
with this kB; however, agreement with experimental points is not so good 
(rVn.d.f. = 4.3). 
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(3) Quenching factors for a particles in CaF 2 (Eu), obtained by the same 
group as in Ref. [36], can be derived from their paper [37]. Alpha peak from 
216 Po (E a = 6.778 MeV [35]) was observed at energy of 1.3 MeV, from 214 Po 
(E a = 7.687 MeV) - at 1.6 MeV, and from 212 Po (E a = 8.784 MeV) - at 
2.0 MeV; thus corresponding quenching factors are equal 0.192, 0.208 and 
0.228, respectively. Because we could expect the same (or similar) condi- 
tions of measurements and data treatment in both [36J and [37], the value 
of kB = 5.1 x 10~ 3 g MeV" 1 cm -2 could fit the data for a's as well. Such a 
curve is presented in Fig. 5c and is in good agreement with the experimental 
data. 



3.2.3 PbW0 4 

Quenching factors for a particles in PbWC>4 crystal scintillator (p = 8.28 g 
cm -3 ) were studied in [38] in the range of energies E a = 2.1 — 5.3 MeV. They 
are shown in Fig. 6 together with fit by Eq. (5) with kB = 10.5 x 10~ 3 g 
MeV" 1 cm" 2 (rVn.d.f. = 3.3). 




a particles 

0.1 - 
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Fig. 6. Quenching factors for a particles in PbWC>4 measured in [38J and their fit 
with kB = 10.5 x 10~ 3 g MeV -1 cm" 2 . 



3.2.4 ZnWO A 

(1) ZnWC>4 (p = 7.8 g cm -3 ) is one of perspective scintillators in searches for 
dark matter particles (see e.g. 6 30 ). Quenching factors for a particles were 
studied in [40p°1. They are shown in Fig. 7a together with fit by Eq. (5) with 
kB = 9.0 x 10~ 3 g MeV -1 cm -2 ; the agreement between the calculated curve 
and the experimental data is good (x 2 /n.d.f. = 0.93). 

10 As in CdW04, quenching factors in ZnWC^ depend on direction of movement 
of a particle; we use here data for direction perpendicular to (010) crystal plane 
measured in wider energy range [40]. 
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Fig. 7. (Color online) (a) Quenching factors for a particles in ZnW04 measured in 
[30] and their fit with the Birks factor kB = 9.0 x 10~ 3 g MeV -1 cm 2 . (b) Quench- 
ing factors for O, Zn and W ions in ZnW04 calculated with this kB value together 
with experimental data for mixture of O, Zn, W factors measured at temperature 
~ 20 mK ©. 

(2) In Fig. 7b, calculated quenching factors for O, Zn and W ions obtained 
with this kB value are presented. Also shown are experimental data collected 
with ZnWC>4 detector in bolometric measurements [6] . Comparing calculations 
with these data, we have to remember that: (a) the latter are in fact data for 
some mixture of quenching factors for O, Zn and W ions; (b) predictions are 
given for kB = 9.0 x lCT 3 g MeV" 1 cm derived by fitting data for a particles 
measured at room temperature while in Ref. [6] temperature was ~ 20 mK. 
To obtain reliable predictions, data for all particles should be collected at the 
same conditions (and with the same data treatment). 



3.2.5 CaWO A 

Quenching factors for CaWC>4 scintillator (p = 6.06 g cm -3 ) probably are the 
most extensively investigated, in particular, because of numerous studies in 
the CRESST experimental searches for dark matter [41]. 

(1) Values for a particles were measured with external and internal a sources 
in [H] in the energy range E a ~ 0.5 — 8 MeV. These data are shown in Fig. 
8 together with their fit by Eq. (5) with kB = 6.2 x 10~ 3 g MeV -1 cm -2 
(x7 n -cl-f. = 6.4). 

(2) Energy dependence of quenching factors for O, Ca and W ions in CaW0 4 
was measured with monoenergetic neutron beam in [42J (for W, only upper 
limits of Qi were obtained, see Fig. 9a). Data for O ions were fitted by Eq. 
(5); the obtained value of the Birks factor is: kB = 8.0 x 10~ 3 g MeV -1 cm -2 
(no surprise that it is different from that determined above for a particles due 
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Fig. 8. Quenching factors for a particles in CaWC>4 measured in |14j and their fit 
with kB = 6.2 x 1(T 3 g MeV" 1 cm" 2 . 

to different conditions of measurements and data treatment I""! Now, fixing 
this kB value, we calculate quenching factors for Ca and W ions; all results 
are shown in Fig. 9a. One can say that calculations for O ions are in reliable 
agreement with the experimental points (x 2 /n.d.f. = 4.8), and also curve for 
W ions is not in contradiction with the measured W limits. 
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Fig. 9. (Color online) Quenching factors for O (circles), Ca (triangles) and W 
(squares) ions in CaW04 measured in [12] and their fits: (a) with the ion total 
SP and kB = 8.0 x 1(T 3 g MeV" 1 cm , (b) with only electronic part of the SP 
(kB = 9.8 x 10" 3 g MeV" 1 cm" 2 ). 

(3) We want to return here to discussion on which stopping power (SP) - total 

11 Formally, x 2 /n.d.f. value for O ions is very high (28.3) that is related with the 
small error bars in the experimental data. However, it should be noted that devia- 
tions of the calculated values from the experimental ones are in the range of only 2 
- 12%. 
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or only electronic part - is better to use in fitting experimental data. Using the 
same ideology, the experimental results for O ions were fitted by Eq. (5) where 
only electronic part of ion SP was taken into account. Obtained Birks factor 
was equal kB = 9.8 x 10 -3 g MeV -1 cm -2 ; fit is worse than that with the total 
SP (see Fig. 9b). With this kB value, quenching curves for Ca and W ions 
were calculated; they are in evident disagreement with the experimental data. 
Thus, use of the ion total SP in fitting and predicting quenching factors in 
the proposed approach allows to describe experimental data in a much better 
way. 

(4) Dependence of the light output and quenching factors in CaWC>4 on atomic 
mass of ion was measured in Ref. [43] by impinging various ions - from H to 
Au - onto the scintillating crystal. All ions had energy of 18 keV. Inverse 
values to the relative light outputs of Eq. (6), normalized to that of electrons 
at 6 keV, were calculated and presented in Table 4 of Ref. 



R' = 1 _ L e(E )/E 
1 Ri(E) U{E)/E 1 ] 



with E — 18 keV and Eq = 6 keV; values of R[ changed from ~ 2 for H 
to ~ 40 for Au. Remembering that the kB values could be different under 
different experimental conditions and data treatment, we will use R[ for H 
to determine kB in this particular measurements and will use this kB value 
to calculate R[ for all other ions. Value of kB = 17 x 10~ 3 g MeV -1 cm -2 
gives R' n = 2.17 that well reproduces values for H given in Table 4 of [43] 
(2.15 ± 0.02 for collection of the scintillation signal during At = 40 /is, and 
2.18 ± 0.02 for At = 50 /as). R[ values for other ions calculated with this kB 



value are shown in Fig. 10a together with experimental results 12 1 of Ref. [43] 
in dependence on ion's Z number. Theoretical points lay on smooth curve well 
fitted by polynomial R\{Z) = a + bZ + cZ 2 with a = 0.87098, b = 0.98708 and 
c = —5.9896 x 10 -3 (also shown in Fig. 10a). 

Comparison of the calculated and experimental values gives extremely high 
value of xV n -d.f. = 112 that is related mainly with deviations for O, Si, 
Ca, Cu and especially Y points. It should be noted that the value for Y 
evidently deviates from the general trend; also error bars in [43] could be 
underestimated. F.e., new data for various ions are presented in Ref. [44] where 
points for O, Si, Ca, Cu and Y are much lower (and in much better agreement 
with the calculated curve), while points for other ions are approximately as 
in Fig. 10a. However, because the data of Ref. [44] were preliminary and still 
are not explained in detail, we do not use them here. 



12 We use data for At = 50 lis that gives more complete collection of scintillating 
signal; however values for At = 40 lis are close. 
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Fig. 10. (Color online) (a) Dependence of inverse of the relative light output 1/Ri 
at Ei = 18 keV, normalized to that for electron at E e = 6 keV on ion's Z number: 
squares are experimental points from [33], and circles are calculated values with 
kB = 17 x 1(T 3 g MeV" 1 cm . (b) Relative light output for a particles with 
kB = 7.5 x 10~ 3 g MeV -1 cm -2 to fit experimental point at 2.3 MeV measured 
at 7 mK in [41]. (c) Predicted curve for Pb ions with kB the same as in (b) and 
experimental point from [TjQ also measured at 7 mK. 

(5) Fig. 14 of Ref. [43J includes also values for a particles at 2.3 MeV 
( 147 Sm) and for Pb ions at 104 keV (nuclear recoil after a decay of 210 Po) 
measured in [JT]. Taking into account that (a) quenching factors are energy 
dependent and (b) data for a and Pb were measured at a temperature 7 mK 
[4~T] . one could not expect perfect agreement between these R\ and R[ for 
all other ions taken in [43J at 18 keV and measured at a room temperature. 
However, inside our ideology, data for a and Pb, taken at the same conditions, 
should be self-consistent. Once more, we can use the relative light output for 
a particle L n (E)/E at 2.3 MeV normalized to that for electrons at 122 keV 
(as in Ref. [41) 13 1 ) to determine the Birks factor; it gives kB = 7.5 x 10~ 3 g 
MeV -1 cm -2 (see Fig. 10b). Curve with this kB value for Pb ions is shown in 
Fig. 10c; at energy of 104 keV it agrees with quenching factor given in [4Hl4'3"] 
(x 2 /n.d.f. = 0.78). Taking into account big difference in atomic numbers of a 
particle and Pb ion (2 and 82, respectively), as well as in their energies (2.3 
MeV and 104 keV), this example demonstrates consistency in description of 
such diverse data. 



3. 2. 6 CsI(Tl ) and CsI(Na ) 

(1) Experimental data for quenching factors of Cs and I ions in CsI(Tl) crystal 
scintillators (p = 4.51 g cm -3 ) published in [45], as well as in [16] (KIMS 



However, numerically it is very close to quenching factor Qi, see Eqs. (5-7). 
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collaboration^ and [47J (TEXONO collaboration) are shown in Fig. 11a, 
lib and 11c, respectively; total energy range of ions was 7 to 135 keV. All 
these data sets are well described by Eq. (5) with the same value of the Birks 
factor: kB = 3.2 x 10" 3 g MeV" 1 cm" 2 (x 2 /n.d.f. value is 1-9, °- 51 and °- 49 , 
respectively). Cs and I nuclei have very close atomic numbers and masses, and 
their quenching factors also are very close (see Fig. 11a). 
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Fig. 11. (Color online) Quenching factors (squares) for Cs and I ions in CsI(Tl) mea- 
sured in: [55] (a), [36] (b), [U] (c) and their fit with the Birks factor kB = 3.2 x 1CT 3 
g MeV" 1 cm -2 ; (d) quenching factors for a particles in CsI(Tl) measured in [58] 
(triangles) and [49] (squares) and calculated quenching curves with kB = 3.2 x 10" 3 
g MeV" 1 cm" 2 and kB = 2.3 x 10" 3 g MeV -1 cm" 2 . 

(2) The KIMS and TEXONO collaborations measured also quenching fac- 
tors for a particles derived from studies of internal contamination of CsI(Tl) 
detectors [4"8ll4"9] . Their values are consistent; see Fig. lid where they are pre- 



sented together with a calculated quenching curve for a particles with the 
same kB = 3.2 x 10" 3 g MeV -1 cm -2 as above for Cs and I ions. However, 
this curve is lower than the experimental points which are much better de- 
scribed by curve with kB = 2.3 x 10" 3 g MeV -1 cm" 2 . This could be some 
underestimation of Q a values in our approach but as well we have to remember 
that the data for Cs and I were taken in devoted measurements with neutron 
beams while data for a particles were collected in separate measurements with 
probably non-identical experimental conditions and details of data treatment 
(temperature in the Cs/I measurements by KIMS was 24.5°C in 
26 — 29°C in a measurements |48j . etc.). 



4"6l. and 



(3) Experimental quenching factors for Cs and I ions in CsI(Na) crystal scin- 
tillators measured in [47)] are shown in Fig. 12a (data for 0.0188% of Na were 
taken but results for other Na amounts are similar). Data points are fitted 



14 In Ref. [46] . data with Tl concentration of 0.128% were chosen but values with 
other Tl concentrations are similar. 
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by Eq. (5) calculated with kB = 5.5 x 1(T 3 g MeV" 1 cm" 2 (x 2 /n.d.f. = 1.9). 
Quenching factors for a particles obtained with this kB value are shown in 
Fig. 12b. 
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Fig. 12. (Color online) (a) Quenching factors for Cs and I ions in CsI(Na) measured 
in [36] and their fit by Eq. (5) with kB = 5.5 x 1CT 3 g MeV -1 cm -2 , (b) Quenching 
factors for a particles in CsI(Na) calculated with the kB value as in (a). 



3.2.7 Nal(Tl) 

(1) Quenching factors for Na and I nuclei in Nal(Tl) crystal scintillators 
(p = 3.67 g cm -3 ) measured in Ref. [50] in energy range of 7 — 215 keV and 
13 — 54 keV, respectively, are presented in Fig. 13a. Value of the Birks factor 
which allows to describe curve for Na ions with Eq. (5) is equal: kB = 3.8 x 10~ 3 
g MeV -1 cm -2 (x 2 /n.d.f. = 1.2). Once found, it also allows to calculate quench- 
ing curve for I ions, and this curve is in excellent agreement (x 2 /n.d.f. = 0.26) 
with the experimental data, as one can see in Fig. 13a. 

(2) In recent work [8], quenching factors for Na ions in Nal(Tl) were measured 
in more detail in low energy region (10 — 100 keV). Obtained experimental 
values are lower than those in [50] . This demands higher kB factor to describe 
higher quenching: kB = 6.5 x 10 -3 g MeV -1 cm~ 2 . Corresponding calculated 
curve well describes experimental values (except of the last point at 100 keV), 
see Fig. 13b (% 2 /n.d.f. is equal 4.0 for all the points and 0.56 with the last 
point excluded). Calculations for I ions are also presented (however quenching 
factors for I were not measured in [8]). 

It is worth noting that with the kB values as in Figs. 13a and 13b, our sug- 
gested procedure predicts a quenching factor for 6 MeV a particles of 0.35 and 
0.25, respectively; these values are lower by a factor 1.5-2 than those typi- 
cally measured (however, under different experimental conditions) for Nal(Tl) 
(see f.e. [EE]). 
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Fig. 13. (Color online) (a) Quenching factors for Na (squares) and I (triangles) in 
Nal(Tl) measured in [50] and calculated curves with kB = 3.8 x 1(T 3 g MeV" 1 
cm -2 , (b) Qi values for Na ions from [8] and fitting curve with kB = 6.5 x 10~ 3 g 
MeV -1 cm~ 2 . Prediction for I ions is also shown. 

(3) Quenching factors for Na and I ions in Nal(Tl) scintillators were also 
measured in other works with monoenergetic neutron beams as: 

- Q Na (£ ^ 18 - 800 keV) ~ 0.30, and Qi(E ~ 20 - 120 keV) ~ 0.1 (with 
some increase at lower energies) [52] ; 

-Q^(E ~ 18-74 keV) = 0.25±0.03, and Q^E ~ 40-100 keV) = 0.08±0.02 

[53]; 

- Qn&(E ~ 50 — 336 keV), values are consistent with constant of 0.27 ± 0.02 

M- 

In addition, in measurements with 252 Cf source quenching factors in some 
effective energy range were obtained as: 

- <2Na(£ = 5-100 keV) = 0.40±0.20, and Qi(E ~ 40-300 keV) = 0.05±0.02 
[5S]; 

- Qn^E = 7-100 keV) = 0.30, and Qi(E ~ 20 - 330 keV) = 0.09 [SB]. 

Comparing energy ranges investigated in these works with behaviour of quench- 
ing curves in Figs. 13a,b, one can note that the energy thresholds were not low 
enough to observe increase of quenching factors at lower energies which is pre- 
dicted in our approach, and at higher energies Qi are consistent with constant 
values (taking into account experimental uncertainties). However, works [81150] 
with lower thresholds give experimental values consistent with this prediction. 

(4) In Ref . [ST] , energy dependence of quenching factor for a particles in range 
of energies E a = 5.7 — 6.8 MeV (from internal contamination of one of Nal(Tl) 
crystals) was obtained as: Q a (E a ) = 0.467(6) + 0.0257(10) x E a , where E a is 
in MeV. This energy dependence (shown in Fig. 14a) can be reproduced by 
Eq. (5) with the Birks factor kB = 1.25 x 10" 3 g MeV" 1 cm 2 much lower 
than those found for data of Refs. [8|5D] in Figs. 13a,b. 
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Fig. 14. (Color online) (a) Quenching factor for a particles measured in [57] (filled 
area represents corridor of uncertainties) and its fit with kB = 1.25 x 10~ 3 g MeV -1 
cm -2 , (b) Predicted quenching curves for Na and I ions with kB = 1.25 x 10~ 3 g 
MeV^ 1 cm~ 2 . Note that the calculated QNa and Q\ are ~ 2 higher than those 
measured in the DAMA and other experiments - however, under different con- 
ditions. 

Following our method, we predict quenching factors for Na and I ions in this 
crystal with the obtained kB value; they are shown in Fig. 14b and are sig- 
nificantly higher (around 0.6 for Na, and 0.3 for I) than the ones measured by 
the DAMA group in [5B] (0.30 and 0.09, respectively), which are also similar 
to other determinations available in literature. 



Considering that our suggested procedure predicts: i) for the cases of Refs. 
[8|50] a quenching factor for a particles in the MeV region lower by a factor 
1.5-2 than those typically measured; ii) quenching factors for Na and I ions 
at low energy higher by a factor about 2 than those measured by |56j and 
other experiments; one could consider the obtained here results with some 
cautious attitude. However, we want to once more remind that quenching fac- 
tors for Na/I ions and a particles were measured in both cases under different 
experimental conditions. It would be of great interest to measure them simul- 
taneously, with the same data taking and treatment. Fig. 10 for CaW04 and 
Fig. 13a for Nal(Tl) give examples that, if such a condition is fulfilled, the 
method gives self-consistent description of data for ions with very different 
(A, Z) values. 

As general consideration, it is worth noting that higher quenching factors 
would always be of big importance for searches of dark matter particles be- 
cause, with a fixed energy threshold of a detector, higher quenching factors 
allow to study lower energies of DM particles. For example, when assuming 
some cases for WIMPs, the energy distribution would drop quasi-exponent ially 
with energy and, thus, this would lead to higher experimental sensitivities in 
the DM searches. 
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3.2.8 CeF 3 



Quenching factors for a particles in CeF 3 crystal scintillator (p = 6.16 g cm -3 ) 
were studied in [58J in the range of energies E a = 2.1 — 8.8 MeV. They are 
shown in Fig. 15 together with fit by Eq. (5) with kB = 11.1 x 10" 3 g MeV -1 
cm" 2 (x 2 /n.d.f. = 1.8). 
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Fig. 15. (Color online) Quenching factors for a particles in CeF3 measured in [58] 
and their fit with kB = 11.1 x 10" 3 g MeV" 1 cm" 2 . 

3.3 Liquid noble gases 
3.3.1 Liquid Xe 

(1) Some data on quenching factors in liquid Xe scintillator (LXe, p = 3.52 
g cm" 3 at —109° C [59j~^l) also support predicted in the described approach 
increase of Qi values at low energies. Results obtained in Ref. [SD] are shown 
in Fig. 16a; they are well described by Eq. (5) with the Birks factor kB = 
3.5 x 10 -4 g MeV -1 cm -2 (x 2 /n.d.f. = 0.36). On the contrary, the Lindhard's 
theory [18] predicts decrease of Qi values at low energies; it is also shown in 
Fig. 16a (calculated using description in [20]). 

(2) Quenching factors for Xe ions in LXe measured in are shown in 
Fig. 16b. They are much lower than those of Ref. |60j, and to reproduce these 



5 Density of liquid Xe depends on pressure and temperature. However, numerical 
value of the Birks factor kB does not depend on value of density if kB is measured 
in g MeV" 1 cm -2 . 

16 More exactly, relative light outputs Li(E)/E normalized to that for electrons 
of 122 keV were measured in [61]; however they are very close here to quenching 
factors, see Eqs. (5-7). 
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Fig. 16. (Color online) Quenching factors for Xe ions in LXe: (a) measured in [UD] 
and fitting curve with kB = 3.5 x 10~ 4 g MeV -1 cm~ 2 ; (b) data from [61] and fitting 

(c) measured in [62] and calculated 



curve with kB = 1.7 x lCM g MeV cm , 
curve with kB = 5.5 x 10 -4 g MeV -1 cm -2 . Lindhard's theoretical prediction is 
shown by dashed line in (a)-(c). (d) Q a quenching curves calculated with different 
kB values together with experimental point for 5.3 MeV a particle from [63J. 



data with high quenching, the Birks factor also should have bigger value: 
kB = 1.7 x lCr 3 g MeV" 1 cm" 2 (x 2 /n.d.f. = 5.0). 

(3) Experimental data from |62j, measured in not so wide energy range (35—69 
keV), are consistent with a constant value (Fig. 16c). However, as well they 
are in good agreement with description by Eq. (5) with the Birks factor kB = 
5.5 x lCT 4 g MeV" 1 cm" 2 (x 2 /n.d.f. = 0.15). 

(4) Quenching curves for a particles calculated with the three above quoted kB 
values are shown in Fig. 16d together with experimental result of Ref. |63j for 
5.3 MeV a particles: Q a = 0.78 ±0.08. Because data of Ref. [63] were obtained 
in experimental conditions different from those in [60|61f62j . no surprise that 
not all curves are in good agreement with the experimental point: to be in 
agreement, quenching factors for all particles should be obtained with the 
same detector and in the same experimental conditions and data treatment. 

(5) It should be noted also that there are other experimental data sets for 
quenching factors of Xe ions in LXe which, on contrary, demonstrate decrease 
of Qi values at lower energies [64]. Results of Refs. [65166] were measured in 
energy range of ~ 47 — 110 keV and ~ 43 — 65 keV, respectively, and are not 
far from constant values, as also could be expected from Figs. 16a-c. 
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3.3.2 Liquid Ar 



Data for liquid argon (p = 1.40 g cm -3 ) are rather scarce. In fact, only one 
measurement was performed in the framework of the WARP project [67] where 
quenching factor for Ar ions was obtained with monoenergetic neutrons of 14 
MeV, however without fixing the angle of neutron scattering. Thus quenching 
factor was obtained for some effective energy range, with mean value of 65 keV, 
as: Q Ar (E ~ 65 keV) = 0.28 ±0.03 (at applied electric field of 1 kV/cm). This 
point can be reproduced by Eq. (5) with the Birks factor of kB = 1.25 x 10~ 3 



g MeV cm" 
theory [T81I20] 



; see Fig. 
is shown. 



17a, where also curve calculated in the Lindhard's 
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Fig. 17. (Color online) (a) Quenching factor for Ar ions in LAr measured in [6 
and fitting curve with kB = 1.25 x 10~ 3 g MeV" 1 cm . Lindhard's theoretical 
prediction is shown by dashed line, (b) Q a quenching curve for LAr calculated with 
kB = 1.25 x 10~ 3 g MeV" 1 cm and experimental data measured in [68] (triangles) 
and [69] (squares). 



Dependence of quenching factor for a particles on energy, calculated with 
this kB value, is shown in Fig. 17b. Experimental values of Q a measured 
in experiments [68.69J also are given; however, because they were obtained in 
different experimental conditions, they are not obliged to lay on the calculated 
curve with kB = 1.25 x 10~ 3 g MeV -1 cm -2 , and even are not obliged to be 
in agreement between themselves (as one can see in Fig. 17b). 



3.3.3 Liquid Ne 

Experimental situation for liquid neon (p = 1.21 g cm -3 ) is even worse than 
for LAr: to-date is only one experiment [11] where the relative light output 
normalized to that for electrons of Eq = 511 keV was measured at the energy 
E = 387 ± 11 keV as: (L Nc (E)/E)/(L e (E )/E ) = 0.26 ± 0.03. Calculated 
quenching curve of Eq. (6) which is normalized to this experimental point is 
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shown in Fig. 18a (kB = 2.0 x 10 3 g MeV 1 cm 2 ); Fig. 18b shows prediction 
for quenching for a particles calculated with this kB factor. 
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Fig. 18. (Color online) (a) Relative light output for Ne ions in liquid Ne. Curve is 
normalized to experimental value at 387 keV from Ref. [11] (kB = 2.0 x 10" 3 g 
MeV -1 cm -2 ), (b) Predicted quenching curve for a's in LNe with kB = 2.0 x 10" 3 
g MeV" 1 cm" 2 . 



3.4 Estimation of quenching factors for nuclear recoils in some scintillators 



While quenching factors for nuclear recoils that constitute such scintillators as 
C 9 H 12 , CaF 2 (Eu), ZnW0 4 , CaW0 4 , CsI(Tl), Nal(Tl), LXe, LAr, LNe can be 
found above, for some other scintillators, which are considered as perspective 
detectors in the DM searches, quenching factors were not measured to-date. 
Below we give estimation of Qi values for nuclear recoils in CdW0 4 , PbW0 4 , 
CeF 3 , Bi 4 Ge 3 0i 2 , LiF and ZnSe scintillators which are based on measured 
quenching factors for a particles in these materials. We should remember, 
of course, that the kB values and thus quenching factors for a particles and 
recoils can be different for the same material in different conditions of measure- 
ments and data treatment. Nevertheless, results given below could be useful 
as providing some initial values of Qi. 

(1) Quenching factors for O, Cd and W ions in CdW0 4 scintillator are shown 
in Fig. 19a for two extreme kB values: kB = 10.1 x 10" 3 g MeV -1 cm" 2 and 
kB = 21.5 x 10" 3 g MeV" 1 cm" 2 . The first kB value was obtained by fitting 
experimental data of [13] for a particles (see Fig. 4a), and the second one by 
fitting data of [34] for a particles and protons (Fig. 4b). 

(2) To calculate Qi for O, W and Pb ions in PbW0 4 (Fig. 19b), we use the 
value of kB = 10.5 x 10" 3 g MeV" 1 cm" 2 obtained by fitting data of [38] for 
a particles (see Fig. 6). 
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Fig. 19. (Color online) Quenching factors for nuclear recoils in scintillators: (a) O, 
Cd and W ions in CdWO/i; (b) O, W and Pb ions in PbWCU; (c) F and Ce ions in 
CeF3; (d) O, Ge and Bi ions in Bi4Ge30i2; (e) Li and F ions in LiF; (f) Zn and Se 
ions in ZnSe. 
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(3) Quenching factors for F and Ce nuclear recoils in CeF3 scintillator (Fig. 
19c) are calculated supposing kB = 11.1 x 10~ 3 g MeV -1 cm -2 obtained by 
fitting data of [5S] for a particles (see Fig. 15). 

(4) For Bi 4 Ge30i2 crystal scintillator (p = 7.13 g cm -3 ), not only quenching 
factors for nuclear recoils were not measured, but also quenching for a particles 
was not studied in detail. Q a values were quoted in several works, but mainly 
for ~ 5.5 MeV a particles from 241 Am source. The results for E a = 5.5 MeV 
are quite different: Q a = 0.17 at 20 mK temperature [70], Q a = 0.20 — 0.21 at 
a room temperature J2B], and Q a = 0.45 at 12 mK temperature [71]. Result 
of Q a = 0.30 ± 0.03 for E a = 5.5 - 8.8 MeV [72J is also known. 

We give in Fig. 19d calculations of quenching factors for O, Ge and Bi recoils 
in Bi4GesOi2 scintillator for two extreme kB values: kB = 2.7 x 10~ 3 g MeV -1 
cm -2 (which reproduces value of Q a = 0.45 for 5.5 MeV a particle of [71] ) 
and kB = 10.9 x 10~ 3 g MeV -1 cm -2 (which reproduces Q a = 0.17 for 5.5 
MeV a particle [70]). 

(5) Information on quenching in LiF crystal scintillator (p = 2.64 g cm -3 ) is 
extremely scarce. We were able to find only one paper where value of Q a = 0.29 
was measured for 5.5 MeV a particles [70 . This value can be reproduced by 
Eq. (5) with the Birks factor kB = 2.5 x 10~ 3 g MeV -1 cm -2 . Quenching 
factors for Li and F ions in LiF with this kB are shown in Fig. 19e. 

(6) ZnSe crystal scintillator (pure and doped by various elements; p = 5.65 
g cm -3 ) is very interesting material in which an extremely low quenching is 
observed for a particles: measured values of Q a are close to 1. For example, 
Q a = 1.0 ± 0.1 for 5.5 MeV a particles from 241 Am was obtained in [73J. 

Several samples of ZnSe doped by O, Al, Cd, Te were investigated in Ref. [74] . 
Q a values for a particles of ~ 5.2 MeV ( 239 Pu) were different for crystals with 
different dopants: for time of a signal collection At = 12.8 ps, value of Q a = 
0.70 was obtained for ZnSe(Cd), and Q a = 0.82 was obtained for ZnSe(O). For 
ZnSe(Te), which is a slow scintillator (with decay time for different components 
of scintillating signal as 30 — 80 ps [74]). even values of Q a > 1 were obtained 
(Q a = 1-25 with At = 0.6 ps, and 1.13 with At = 12.8 ps). Evidently even 
collection time of At = 12.8 ps is not long enough to collect total signal in 
this slow scintillator. We could suppose that this is a reason of values Q a > 1, 
and collection of a signal during longer times would give an usual situation 
with Q a < S 

At bolometric temperatures, when signals in ZnSe are extremely long (> 150 



17 Similar situation with protons which have faster signals than 7 quanta in CsI(Tl) 
was already mentioned in section 2: Q p > 1 was obtained for At = 1 ps, and Q p < 1 
for At = 7 ps. 
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ms), it was possible even to obtain values of Q a ~ 4 |75j. 

We have to note here that Eq. (5) supposes quenching of a signal from highly 
ionizing particle in comparison with that from electrons, and it is impossible 
to describe any enhancement on its basis (staying with physical values of 
kB > 0). For kB = 0, Eq. (5) gives quenching factor Qi = 1 for any particle. 
In Fig. 19f, we suppose usual situation when a signal is totally collected and 
Qi < 1. Quenching factors are drawn for Zn and Se ions in ZnSe with kB 
value as: kB = 1.9 x KT 5 g MeV" 1 cm" 2 (which gives Q a = 0.99 for 5.5 MeV 
a particles) and kB = 8.4 x 1(T 4 g MeV" 1 cm" 2 (Q a = 0.70 for 5.5 MeV a 
particles) . 



4 Conclusions 



Semi-empirical and quite simple in realization method of calculation of quench- 
ing factors for scintillators was described in this work. It is based on the clas- 
sical Birks formula with the total stopping powers for electrons and ions, and 
has only one parameter: the Birks factor kB. Value of this factor for a given 
scintillating material can be different in different conditions of measurements 
and data treatment. However, if experimental conditions and treatment of 
data are fixed, hypothesis that kB has the same value for particles of different 
kinds gives reliable results. Once the kB is found by fitting quenching factors 
for particles of one kind and in some range of energies (e.g. for a particles from 
internal contamination of a detector by U/Th chains and/or by 147 Sm, 190 Pt 
with energies of a few MeV), it can be used to calculate quenching factors 
for particles of another kinds and for another energies of interest (e.g. for low 
energy nuclear recoils). Many examples were given for materials which, fur- 
thermore, have different mechanisms of scintillation: organic scintillators (solid 
CsHg, and liquid Chilis, C9H12); crystal scintillators (pure CdWOzi, PbWC>4, 
ZnW0 4 , CaW0 4 , CeF 3 , and doped CaF 2 (Eu), CsI(Tl), CsI(Na), Nal(Tl)); 
and liquid noble gases (LXe). It was demonstrated for many cases that the 
method allows not only to describe measured data for ions of one kind in a re- 
liable way but also to predict behaviour of quenching factors for other particles 
which sometimes is immediately confirmed by already existing experimental 
data - sometimes worse, sometimes better, and sometimes very good, but at 
least in a rough agreement. Some predictions (e.g. for LNe, LiF and others) 
could be checked in near future. 

Stopping powers for electrons and ions are calculated with the ESTAR and 
SRIM codes, respectively, which in fact present to-date state-of-art software 
in this field. It is easy to use these programs and they are publicly available; 
this makes Qi calculations quite simple. 
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Calculations with the SRIM package have some tendency to overestimate 
quenching factors for a particles at energies around ~ 2 MeV and underesti- 
mate them at high energies (> 8 MeV) as can be seen in Fig. 4a for CdW0 4 , 
Fig. 5a for CaF 2 (Eu), Fig. 8 for CaW0 4 , and Fig. 15 for CeF 3 . At the same 
time, calculation of the stopping powers for a particles with the ASTAR pack- 
age gave better description of Q a in CaF2(Eu) scintillator at lower energies 
(see Fig. 5a). For some other materials difference between ASTAR and SRIM 
calculations was not big (see Fig. 2 for CsI(Tl) and Fig. 3a for C 8 H 8 ). Evi- 
dently Qi values will depend on how one calculates stopping powers for ions 
and electrons, and it is a pity that stopping powers could not be computed in 
framework of the same package for any particle (SRIM calculates SP for ions 
in any substance but does not calculate SP for electrons; and STAR gives SP 
for electrons in any material but SP for ions are possible only for protons and 
a particles and for a limited list of materials). 

Quenching factors calculated in the presented approach in general increase 
at low energies, and this encourages experimental searches for dark mat- 
ter particles. Estimations of quenching factors for nuclear recoils are given 
for some scintillators where experimental data are absent (CdW0 4 , PbW0 4 , 
CeF 3 , Bi 4 Ge 3 0i 2 , LiF, ZnSe). 
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